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Abstract 

We complete the one loop renormalization of the HQET lagrangian at 
©(l/mg) including four fermion operators with two heavy and two light 
quark fields in the operator basis. It is shown that as a consequence the 
short distance coefficients of the operators bilinear in the heavy quark 
field receive nontrivial corrections. 
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1 Introduction 



Heavy Quark Effective Field Theory (HQET) [Q has become a weU established 
theoretical tool for the description of hadrons containing one heavy quark Q . 
This derives from the fact that it is a systematic expansion in inverse powers of 
the heavy quark mass toq with well defined and calculable coefficients. Further- 
more, its realization of the spin and flavor symmetry of the low energy theory 
is a phenomenologically powerful tool. The l/mg expansion has already been 
applied sucessfuUy to phenomenological problems such as the determination of 

Besides its phenomenological application HQET posseses interesting the- 
oretical features which already show up if one studies the HQET lagrangian 
itself. On one hand it is relevant for power corrections to the 1/toq expansion 
of hadron masses and on the other hand allows for a study of reparametrization 
invariance [ p^ . 

However, the complete construction of HQET has to include radiative correc- 
tions. These are computable order by order in perturbation theory by matching 
HQET to full QCD at a perturbative scale /i = mq followed by the renormal- 
ization group running of the respective operators. As indespensable ingredients 
both the coefficients at the matching scale and the anomalous dimensions of the 
operators are needed. 

In this paper we concentrate on the one loop renormalization of the la- 
grangian at 0{l/mQ). Results for the terms of 0{l/mQ) and ©(l/mg) as 
well as a matching calculation for the terms of 0(1/171%) are already known 

In |1S| the one loop anomalous dimensions of the effective lagrangian at 
0(l/mQ) have been calculated. However, these results were incomplete in the 
sense that only operators bilinear in the heavy quark field had been included 
in the operator basis. In the present paper we complete the renormalization of 
the effective lagrangian extending the operator basis to four fermion operators 
composed of two heavy and two light fields. At ©(l/mg) it is well known that 
the running of the short-distance coefficient of the Darwin operator is modified 
in the presence of such operators. At higher orders this effect should be expected 
to cause more drastical consequences because of the large operator bases. In 
fact, it will be shown that at 0(1/771^) the coefficients of the bilinear operators 
receive nontrivial corrections in the presence of heavy-light operators. 

A possible phenomenological application of our results is the estimate of 
0(l/mQ)-corrections in the operator product expansion of the B-meson decay 
width. It is well known that to this order the width is significantly affected 
by spectator quark effects induced by heavy-light operators appearing in the 
OPE. However, there are also time ordered products of the effective lagrangian 
with the local operators in the OPE. At this stage heavy-light operators in the 
lagrangian and their short distance corrections which are the subject of this 
paper induce additional spectator effects. 

This note is organized as follows: in section ^we introduce our operator basis 
and discuss its reduction to a set of linearly independent operators in section ^ 
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Our result for the anomalous dimensions is presented in section Finally, we 
present the logarithmic contributions to the short distance coefficients of the 
effective lagrangian in section p| and conclusions in section pi 



2 Operator basis 

The lagrangian of HQET can be written as 

£ = K{ivD)h, + (^^^^^^"^ + '^"sht , (1) 

where 

Aight = -if^.F"^'' + qi^P - m,)q (2) 

describes the dynamics of the light degrees of freedom. For simplicity we con- 
sider only one light quark with mass niq. Additional light flavor are taken care 
of by summing over the index q. 

To lowest order in the heavy mass expansion only one operator, hy{ivD)hy, 
shows up which has the celebrated spin- and flavor-symmety properties. 

However there are power corrections of ©(I/ttiq) to the heavy-quark limit 
which break these symmetries. They are given by a sum of operators of appro- 
priate canonical dimension multiplied by short distance coefficients: 

= (3) 

i 

At given order 0(l/mQ) all operators of the respective canonical dimension 
allowed by the symmetries of the effective theory may contribute. 

First of all there are the operators bilinear in the heavy quark fields already 
appearing in the tree level lagrangian (denoted H-operators) . 

Operators containing two heavy and an even number of additional light 
quark fields appear if one matches QCD amplitudes with heavy and light exter- 
nal fields onto HQET to one loop. Their inclusion accounts for spectator quark 
effects inside hadronic states, between which these operators are sandwiched 
when calculating physical matrix elements. However, since their massdimension 
must be larger or equal to 6 these operators show up at 0{l/mQ) for the first 
time. To this order they can be constructed by two heavy and two light quark 
fields. At 0{l/mQ) a covariant derivative or a factor iriq increases the mass 
dimension. In what follows these operators are referred to as HL-operators. 

Since we remain in the one particle sector of HQET operators with four or 
more heavy quark fields are ommited. 

To complete the operator basis we should also consider operators constructed 
of four light quark fields and the corresponding penguin operator. However, 
to one loop order these operators only mix with the HL-operators. Since we 
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are focusing on the corrections to the H-operators we therefore neglect such 
operators. 

Before we proceed to 0{l/m^) let us recall the operators appearing at lower 
orders. 

At ©(l/mg) we choose the conventional basis 



of) = KiiDfhy 
of) = K{ivDfhy , 



and at 0(l/m^) 



of) = KiDf,{ivD)iD''hy 
of) = K{ivD){iDfhy 



O 



(2) 



hv{ivD)^hv 



of) = Ki(7''^iD^{ivD)iDxhy 

of) = hy{iDf{ivD)K 

of) = hy{ivD)ia^'^iD^iDxK 



(4) 



(5) 



of) = Kia^'^iD^,iDx{ivD)K 



The definition of the covariant derivative iD = id + ggT'^A"- and field strength 
tensor F^^,T°' = —i/gs[iDi^,iD„] follows usual conventions. 
We choose the HL-operators at 0{l/mQ) as 



^i2)s/o _ „2 



9s[QC:/oQ][hvCt/ohv] 



(6) 



M 



(2)s/o _ 



= g^^lqia^^-C^^^h^ta^X^/ohv 



/o 



where C° = 1 in the color singlett and = T" in the color octett case. 

At 0(l/mg) we find 13 local H-operators contributing to physical matrix- 
elements: 



of) 


= hiD^{ivDfiD''h 


of) 


= h{iDf{iDfh 


of) 


= JiiD^iiDfiDi'h 


of) 


= hiD^iDJD^'iD"h 


0(3) 


= hicr^"'iD^,{ivDfiD^h 


0(3) 


= hia'"'iDp,iD4iDfh 


of) 


= Jiia^'^iDpiD^iD^iDPh 


0(3) 


= hiai"'{iDfiDp,iD^h 


of) 


= hi(ji"'iDf,{iDfiD^h 


0(3) 


= hia'"'iD^iDpiDJDPh 


of) 


= Ua^^iDpiD^iDHD^h 


0(3) 


= glhF'^i^^F^.h 


off 


= g^hv^vPF'^^^-'F^ph 







(7) 



Note that, in contrast to the case of the lower order operators, we omit operators 
vanishing by the heavy quark equation of motion (EOM) 



{ivD)hy = 



(8) 
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from the very beginning. This is justified since we are not going to insert such 
operators in time ordered products of order 0(1/toq) or higher. We may not 
remove such operators from the operatorbases at ©(l/mg) and ©(l/mg) since 
time ordered products of these operators at 0{l/mQ) may be identical to local 
physical operators and contribute to their anomalous dimensions (see section 
3.2). However, the following operators are needed as parts of the gluon EOM 
and though being irrelevant we will keep them in the basis: 



O 



(3) 
14 



K{ivDf{iDfh^ 
O'C^ = K{iDf{ivDfK 



)(3) _ 

= h^{ivD)iD^{ivD)iD^'K 



0^5^ = h^{ivD){iDf{ivD)h^ 
Of^ ^hyiDf,{ivD)iD''{ivD)hy 



To construct a basis of HL-operators at 0{l/mQ) we let a covariant deriva- 
tive act either on a light or a heavy quark field. A further classification has to 
be taken according to the direction in which the covariant derivative acts. We 
define 



9sAlT<^ 



(9) 



and it is understood that in the octett case the covariant derivative stands 
left/right of the color matrix when acting to the left /right. With these defini- 
tions the HL-operators of mass dimension 7 are choosen as: 



M 
M 



(3/i)s/o 
2± 

(3/t)s/o 
3± 



(3/t)s/o 
6± 

(3/i)s/o 
7± 



M 
M 



M 



^{3h)s/o 

{3h)s/o 
10± 

{3l)s/o 
4± 

(3i)s/o 
5± 



M 
M 



= ±g^A0C:/MhvC:/,{^vD^)K] 

= ±gl[qta^^\xC:/,q] [h^C^^jD^K] 

= ±g'Aqi5l^C:^oQ] [K^5l^.C:/,{^vD^)h,] 

= ±g',[qj5C:/,q] [Kl5C:/,tp^K] 

= ±g^M:/o^P^<l][hvC:/M 

= ±gl[qC:^,{tvD^)q][hX:/ohv] 

^ ±g'MC:/o{^vD^MhvC:/M 

= ±g^M^^''vxC:/^iDtq] [KC:/,K] 

= ±g^^[qia^''C^/^{ivD^)q] [Kia^,r^/^K] 



(10) 
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_yV^(3/)«/o 

(30s/o 
7± 

{3l)s/o 
8± 



M 
M 



_^(3/)«/o 
. ,(3/)s/o 



±.9?[g75Q<,*i?±g][/i.757''Qo^.] 



To complete the set of local operators we also consider operators of lower 
dimension multiplied by an appropriate power of the light quark mass: 



(1) 



"q^i ' '■'■q^i ' iii-qJ^il (11) 

These operators are needed as counterterms during renormalization as long as 
the light quark is kept massive. However the mass dependent HL-operator 
can be removed with help of the light-quark EOM. This is not the case for 
the mass dependent H-operators which are induced by penguin type diagrams. 
They cause corrections proportional to powers of niq/mQ to the short dis- 
tance coefficients of the respective operators. Note that the operators 
i = 1,2,8,10 vanish by the heavy quark EOM. Instead of disregarding these 
operators from the very beginnig we will keep them in the basis since they show 
up in several operator identities. 

In addition to the local operators there are the time ordered products of 
dimension 7 composed of lower dimension operators: 



(12) 



111) 



i T 



/o(l) ^^(2) 



-r _ 

ijk 

{12hl)s/o _ 



-Sijk T 



iT 



0(l)^_^(2)./o 



1 = 1,2,3 j = l,...,7 
z,j,fc=l,2,3 i<j<k (12) 
1 = 1,2,3 j = l,...,4 



™«(1 - o'^ij)*T 



0(1) 0(1) 



(13) 



The symmetry factor Sijk equals 1, 1/2 or 1/6, if no, two or all inserted operators 
are identical. Again a mass dependent T-product 

1 

2^ 

contributes. Although the short distance coefficients of time ordered products 
are the products of the coefficents of their operator components, local operators 
are required for their renormalization. This in turn modifies the running of the 
coefficents of the local operators. 



3 Reduction of the operator basis 

The operator basis presented in the previous section is overcomplete, i.e. there 
exist several interdependencies between the operators. To arrive at a phys- 
ical basis of linearely independent operators all redundant operators have to 
be removed, otherwise artificial gauge dependencies may show up in the corre- 
sponding short distance coefficients psj. 
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3.1 Momentum conservation 



Since all operators are integrated over their argument, momentum conservation 
is manifest. 

In case of the H-operators the only effect is that the direction in which the 
covariant derivative acts is irrelevant. 

The partial derivative of a HL-operator of mass dimension 6 vanishes. Per- 
forming the differentiation explicitly yields relations among the HL-operators 
of dimension 7: 

= df,[qT''''q][Kr,K] 

= [qiD~r''''q][h,T,K] + [qT^'''iD+q][KT,h,] (14) 
+ [qT^'''q][KlD-T,K] + [qT^'''q][Kr,lD+K] 
= d^[qT,q][hT^''K] 

= [qiD-r,q][h,T''''h,] + [qT,q][hr^''iD+h,] (15) 
+ [qT,q][KtD;^T'^''h,] + [qT,q][Kr>'''iD+K] 

Choosing appropriate Dirac matrices F^j, and this allows us to remove the 
following operators from the basis (in favour of the combination on the r.h.s): 

m[ 

Mi 

Mi+'''" = M):'l'''" -M'f^'f'" + M'f^l!''" (16) 

m\ 

Ml 

(17) 



(3/)s/o 
1+ 






^{3h)s/o 


(3l)s/o 
3+ 


= M^il^'^° - 




^(3h)s/o 


(3/)s/o 
5+ 


= Mfl^'^° - 


Mlf/° + 




:{3l)s/o 
6+ 






^(3h)s/o 


i3l)s/o 
7+ 








(3/)s/o 
10- 


^ -'^'10+ 







3.2 Contraction Identities 

Time ordered products, in which a (w£')ft,^-term acts on an internal line, are 
identical to a combination of local operators: 

iT [hyF{iD){ivD)h^,hyGiiD)hy] = -KF{iD)G{iD)K + . . . (18) 

If only H-operators are involved such contraction identities have been studied 
in |l^ , where it has been shown that all time ordered products of this type 
are redundant and can be removed in favour of local H-operators. This way 
the anomalous dimensions of the local H-operators receive corrections that are 
crucial for the gauge independence of their short distance coefficients. 
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In the presence of HL-operators additional contraction identies arise: 

iT[[KitvD)^K],g^,[qT^^q][hXf,K]\ = -g^A<l^''q][hAivD-)T'^K] 

~gl[qTf^q][KT'^{zvD+)K] (19) 

Since the T-product is only related to HL-operators vanishing by the heavy 
quark EOM the additional contraction identities are irrelevant. 

3.3 Light Quark EOM 

Using the equation of motion for the light quark 

ip^q ~ rriqq — qip^ = niqq (20) 

we can remove the mass dependent HL-operators in (|ll|) from the operator 
basis: 

mg[qr''q][KT'^K] = [qT''ip+ q][hXf,K] (21) 

-[qip-r''q][Kr'^K] 



Taking a symmetrical combination of both equations (|21j) and choosing the dirac 
matrices and T'^ appropriately we derive the relations 

mX'^^^° = iA^f_^)^/°-iA^(f^/°-iAl(3';)^/° + iA^fi)^/° (22) 

,1 ui3h)s/o _ 1 . .(3/)s/o 1 ,.(3l)s/o 

^■'^ho- 2 2 

in which momentum conservation has been already used. 

Subtracting both equations ( |2l] ) yields relations among the operators of di- 
mension 7, written generically as: 

= [qT^'lp+q][hXf,K] + [qip-T^'q][hT'^h,] (23) 

Choosing appropriate dirac matrices and F^ this identity allows us to remove 
the following operators from the basis (in favour of the operators on the r.h.s.): 

M 



M 



, {3h)s/o 








^{3V)s/o 






^i3l)s/o 


, {3h)s/o 

'io+ 




2^(3^)-/° ^ 





7 



3.4 Gluon EOM 

Taking the functional derivative of the lagrangian 

C = -ii^^.F'^'''' + qitp ~ rUg) + K{ivD)K (25) 

with respect to the gluon field A\ yields the gluon EOM 

jjbapat^X ^ g^q^^T'q + g^v^KT^'K = , (26) 

where D"^^ — S'^^'d — gsf'^''^A^ is the covariant derivative in the adjoint represen- 
tation. We multiply T^v\ and sandwich the resulting expression between heavy 
quark spinors to get 

K[iD^, {ivD)]]K = g^AhvT''K][q^q]+g^^[KT''h,][KT''h,] . (27) 

The four-heavy-quark operator on the left hand side can be omitted, since at 
least to one loop order such operators only mix with themselves or appear as 
counterterms and therefore decouple from the renormalization group flow of the 
other operators. 

In our conventions ( p7|) reads 

= -20p^ + of + of . (28) 

To derive an analogous relation among operators of dimension 7 we multiply 
6|) with T'T'^^ and sandwich between —h^iD^ and hy 

- gl[qixT\][KiD-T''T''^K] = ^2KT''^iD,iD^iDxiD''hy 

+KT''^iD^iDx{iDfK (29) 
+KT''^iD^{iDfiDxK 

or between and iD^h^ 

g^,[q-fxT\][hyT'r''\D+K] = -2h,r^''iD^iD,iD^iDxK 

+KT^''{iDfiD^iDxK (30) 
+KT^''iDy{iDfiDxK ■ 

Choosing Ti,x = v^vx, gv\,i<yu\ yields the following relations: 

_^(3/.)o ^ -20'^^ +Of^ +Of^ 
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(n(3) _ 






















-20f^ 







(31) 



Note that the relations involving jWf^^° are consistent with the first relation in 

( p4| ) and the relations involving AA^2±'° with the heavy quark EOM. 

Taking the time ordered product of ( p8|) with a dimension 5 operator 0^-^'' 
yields 

^(i2./)o^_27:(/^)^^(i2)^^(i2)^ (32) 

where the last two T-products on the right hand side obey a contraction identity 
and are related to local operators of'\ Note that in the case 1 = 3 ( p2[ ) is 
consistent with the contraction identity (^9|), if the heavy quark EOM is applied. 

3.5 Reduced Operator basis 

Applying the identities of the previous subsections to the full operator basis 
leaves us with the set of physical operators listed in figure |l|. In terms of this 
set of linearely independent operators the effective lagrangian at ©(l/mg) reads 

£(3) ^CB-B + Co,^-d„, + CM- M + Ct -f . (33) 

The different types of operators and their coefficients are collected in vectors as 
defined in figure |^. 



4 Anomalous dimensions 

In order to compute the anomalous dimensions of the operators (B, Om, ■M^T)^ 
in general one has to calculate the pole parts of all divergent IPI Greensfunctions 
with a certain operator inserted and express the result in terms of tree-level 
operators 

where Ai runs over all operators of the basis. From this one can directly read 
off the anomalous dimensions . 

Since we are working in the background field gauge jl^, |l^, it suffices to 
consider IPI Greensfunctions with cither two external heavy quark legs and, 
depending on the dimension, one or two background fields, or two heavy and 
two light fermion fields. Figure |^ shows the mixing properties of the operator 
basis by diagrammatic examples. The procedure how to take care of the var- 
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operator 


index 


shorthand (number) 


Of 


j = 1, . . . , 13 




ijk 


hj,k^l,2 i<j<k 


B (21) 


^(12) 


= 1,2 






i = 1,2 






i = l,2 


d,n (7) 










t = 3,5,6,7,9 






i = 3,4,5,6,7,9 






i = 5,6,7,9 






i = 5,6,7,9 


M (43) 


^ms 


i = 2,8,9 






1 = 1,... ,9 




_^(30o 


i = 2,8,9 




^(3/)o 
I — 


?: = i,... ,9 




^{12hl)s 
ij 


z = l,2 j = l,...,4 


f (14) 


^{12hl)o 
ij 


z = l,2 J = 1,3,4 



Figure 1: Physical operator basis at O (1/171%). 



ious operator identities consistently follows the lines of |l^. We present our 
result for the anomalous dimension matrix corresponding to the operator basis 
{B, Om,-M,T) as a block matrix: 








IBM 


\ 





7o„o,„ 


lO^M 





1MB 




IMM 





\^1TB 





ITM 


ITT J 



(35) 



The mixing of T with S is a consequence of the gluon EOM removing HL- 
operators in favour of (see (^)), since to one loop order there are no 

penguin diagrams which would require Ol as counterterms. 

Furthermore the anomalous dimensions of the H-operators (i.e. the entry 
Ibb) are affected by the gluon EOM since some of the A^-counterterms may be 
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diagram (example) 



local counterterm 



^(3) 



B 





O. 



(3) 



IBB 





M 



{3h/l)s/o 



IBM 



IBB 





^20(1) 



M 





j^{-ih/l)s/o 



IMM 



T 



j^{3h/l)s/o 



ITM 



M 



{3h/l)s/o _^ .^(3) 



Of' (gluon EOM) 



7Te 



Figure 2: Mixing properties and operator identities. The black blob represents 
an operator of appopriate mass dimension. 

removed using (^l]). This in turn implies corrections to the Wilson coefficients 
of the local operators in B. 

The operators Om require local counterterms of the form mgA^f which 
are removed by means of ( p2[ ) in favour of operators A^^^''^'^''^°, thereby ex- 
plaining ^o,„M- The even more complicated situation that some of the latter 
obey the gluon EOM fortunately does not occure. 

Since the submatrices of the block matrix (|3^) are too large we refrain from 
their presentation. They are available from the author on request. 

5 Renormalization group logarithms 

With the one loop anomalous dimensions and the tree level matching coefficients 
we are now in the position to solve the renormalization group equation for the 
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Wilson coefEcients of the effective lagrangian at ©(I/toq): 

d 



^, C(3)(a.)+7(')Tc(3)(^)^0 (36) 
am II 

Here C^'^^ (/i) denotes collectively the coefficents (Cg, Co„^ , Cm , Ct) of the phys- 
ical operators {B, Om, ■M^'T) in the effective lagrangian (|3^). Since an analytical 
diagonalization of 7*^'^^ seems to be difficult we restrict ourselves to the calcula- 
tion of the first logarithmic correction oc as ln(/i/mQ) in the coefficients C*-^-'(/x). 
The exact solution of (Isl) reads 



with the one loop running coupling 



l_2/3^-^(^^)ln(^) (38) 
as{mQ) TT ruQ 

where /?^°) — (33 — 2n/)/12 in the presence of light flavors. Expanding ( ^7|) 
to first order in the strong coupling we get 

Ct3)(^)=Ct3)(^Q)-(^)ln(^)7(3)Tc(3)(^Q) + 0((^)2). (39) 

TT rUQ TT 

With our result 7(3)^ ^jjj^ ^j^g |;j.gg jgvel matching coefficients C'3)(77jg) the 
Wilson coefficients are easily calculated. Note that the operators in O™, AA 
and T are not present at tree level, i.e. the corresponding entries in C'^'^^tiiq) 
are zero. In figure I the non vanishing coefficients are shown. The remainig 
coefficients are ©(a^) and do not contribute in our approximation. 

The coefficients of the time ordered products are given by the product of the 
coefficients of their operator components. 

(3) 

Note, that the coefiicients CI for i = 6, 8, 9, 10, 11 are modified with respect 
to the value they would get if one disregards the HL-operators (in brackets). 
This correction is a consequence of the gluon EOM which relates certain HL- 
operators to H-operators (see (^)). In general the same is true for the coef- 
ficients C\'^\ i = 2, 3,4 but to one loop order the two operators Ai'f^^° which 
would cause the corrections always appear in the combination A^^^"*" — M'f^^° 



which vanishes identically according to (31) or the first relation in ( 



6 Conclusions 

In this paper we have completed the one loop renormalization of the HQET 
lagrangian at 0{l/rn^). In addition to H-operators which are bilinear in the 
heavy quark field HL-operators consisting of two heavy and two light quark 
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Wilson coefficient 


tree level 


coefficient of (Q!s(/ij/7rj in(/i/TOQj 


„{3) 
'-'1 


2 


-25/3 Ca - 32/3 Cf 


„{3) 
^■2 


-1 


-1/2 Cyi 


„{3) 
'-'3 


-1 


ACa + 8/3 Cf 


„{3) 
O4 


1 




„{3) 
^5 


-2 


5/3 + 8/3Cf 


„{3) 
^6 


1 


1/2 Ca (-Ca) 


„{3) 


1 


-13/6 - 8/3 Cf 


„{3) 


1 


1/2 Ca (-Ca) 


„{3) 


1 


-Ca (-4Ca) 


„{3) 


-1 


3/2 Ca (9/2 Ca) 




-1 


3/2 Ca (9/2 Ca) 


„{3) 
^12 





1/12 Ca 


„(3) 
'-'IS 





-1/3 Ca 


^(3/1) S 

^7+ 





-2 Ca Cf + 1/2 C^ + 2 C| - 1/2 


Cj_ 





-2CaCf + 1/2C| + 2C|. - 1/2 


^{3h)o 
^7+ 





2Ca -4Cf 


Q(3h)o 





Ca - 4 Cf 


(j{3l)s 





-20/3 Ca Cf + 5/3 C^ + 20/3 C| - 5/3 


(j{3l)s 





32/3 Ca Cf - 8/3 C^- 32/3 C| + 8/3 


fj(3l)o 





5 Ca -40/3 Cf 


r<(3i)o 





-8 Ca + 64/3 Cf 


r<{3i)o 

^6- 





-Ca 


q(31)o 





Ca 



Figure 3: Wilson coefficients of tlie physical operator basis. 



fields are included in the operator basis. We have shown that there exist sev- 
eral interdependencies between the operators of the complete operator basis 
which makes its reduction to a basis of linearely independent operators a non- 
trivial task. The anomalous dimensions of the H-operators have already been 
calculated in |l3|. However in the presence of HL-operators the short dis- 
tance coefficents of the H~operators are modified. This is a consequence of the 
gluon EOM which relates certain HL-operators to H-operators thereby mod- 
ifying their anomalous dimensions and coefficients. This effect is well known 
at ©(l/mg) where it causes corrections to the coefficient of the Darwin oper- 
ator 01 . However, this correction is weak of 0{c?g) whereas at ©(I/toq) the 
coefficients receive corrections already at 0(as). 
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